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Abstract 

In this paper we discuss the detection of CP and T-violation effects in the frame¬ 
work of a neutrino factory. We introduce three quantities, which are good dis¬ 
criminants for a non vanishing complex phase (<5) in the 3x3 neutrino mixing 
matrix: A g, A cp and At- We find that these three discriminants (in vacuum) all 
scale with L/E Vl where L is the baseline and E u the neutrino energy. Matter ef¬ 
fects modify the scaling, but these effects are large enough to spoil the sensitivity 
only for baselines larger than 5000 km. So, in the hypothesis of constant neutrino 
factory power (i.e. number of muons inversely proportional to muon energy), the 
sensitivity on the (5-phase is independent of the baseline chosen. Specially inter¬ 
esting is the direct measurement of T-violation from the “wrong-sign” electron 
channel (i.e. the A t discriminant), which involves a comparison of the v e —> ry 
and ry —> v e oscillation rates. However, the ry —> v e measurement requires 
magnetic discrimination of the electron charge, experimentally very challenging 
in a neutrino detector. Since the direction of the electron curvature has to be 
estimated before the start of the electromagnetic shower, low-energy neutrino 
beams and hence short baselines, are preferred. In this paper we show, as an 
example, the exclusion regions in the A m\ 2 ~ <5 plane using the A cp and A t 
discriminants for two concrete cases keeping the same L/E u ratio (730 km/ 7.5 
GeV and 2900 km/30 GeV). We obtain a similar excluded region provided that 
the electron detection efficiency is ~20% and the charge confusion 0.1%. The 
A m\ 2 compatible with the LMA solar data can be tested with a flux of 5xl0 21 
muons. We compare these results with the fit of the visible energy distributions. 
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1 Introduction 


The firmly established disappearance of muon neutrinos of cosmic ray origin [ij strongly 
points toward the existence of neutrino oscillations |J. 

The first generation long baseline (LBL) experiments — K2K [8], MINOS [4|, 
OPERA 0 and ICARUS H [7], || — will search for a conclusive and unambiguous 
signature of the oscillation mechanism using artificial neutrino beams produced by the 
“traditional” meson-decay method. They will provide the first precise measurements 
of the parameters governing the main muon disappearance mechanism. 

In contrast, a neutrino factory |9|, [L4|] is understood as a machine where low energy 
muons of a given charge are accelerated in a storage ring. Neutrino factories raised the 
interest of the physics community, since they appear natural follow-ups to the current 
experimental LBL program and could open the way to future muon colliders. 


The real physics potential of a neutrino factory comes from its ability to test in a 
very clean and high statistics environment all possible flavor oscillation transitions. 
As many studies have shown 0, the physics potential of such facilities are very 
vast. An entry-level neutrino factory could test the LSND signal in a background 
free environment [|TTJ . A neutrino factory source would be of sufficiently high intensity 
to perform very long baseline (transcontinental) experiments. 

The ability to measure all possible neutrino oscillations will provide very stringent 
information on all the elements of the neutrino mixing matrix and on the mass pattern 
of the neutrinos. In a 3 x 3-mixing scenario, the mixing matrix, which should be unitary, 
is determined by three angles and a complex phase. As studies have shown JT0[, precise 


determination of two angles and of the largest mass difference will be obtained. In 
addition, a test of the unitarity of the matrix could be performed. 


Apart from being able to measure very precisely all the magnitude of the elements 
of the mixing matrix, the more challenging and most interesting goal of the neutrino 
factory will be the search for effects related to the complex phase of the mixing matrix. 
The complex phase will in general alter the neutrino flavor oscillation probabilities, and 
will most strikingly introduce a difference of transition probabilities between neutrinos 
and antineutrinos (so called CP-violation effects), and between time-reversed transi¬ 
tions (so called T -violation effects) [^Uj] . Neutrino factories should provide the intense 
and well controlled beams needed to perform these studies. 


Two fundamental experimental parameters for these studies are (1) the energy E 
of the stored muon beam which determines the neutrino energy E v spectrum and (2) 
the baseline L between the source and the detector. In earlier works, one has advocated 
high energy neutrinos, with a consensus around ~ 30 — 50 GeV coupled to long 
baselines L ~ 3000 km. In Ref. |R|, a much lower muon energy in the range P M ~ 1 GeV 
coupled to short baselines L ~ 100 km has been proposed. This apparent contradiction 


was recently addressed by Lipari|16 
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In this paper, we expand our early workJTTj] on the detection of CP-violation at 
a neutrino factory. In particular, we concentrate on the general strategies to detect 
effects related to the complex phase of the neutrino mixing matrix. The complete 
and comprehensive detection of the effect is very difficult for terrestrial experiments, 
as it would require L/E v values simultaneously in the range of solar and atmospheric 
neutrinos. We concentrate on medium L/E u and demonstrate that the relevant effects 
scales as L/E u and that there is an optimal choice for the ratio L/E v , given by the 
neutrino mass pattern. Hence this yields an priori large freedom for the choice of L and 
E v . We show however that for too large L, matter effects destroy the sensitivity to the 
complex phase, so baselines smaller than L ~ 5000 km are preferred. In addition, owing 
to the linear rise of the neutrino cross-section with neutrino energy at high energy, the 
statistical significance of the effect scales with E v , and hence grows linearly with L, for 
L/E u constant. 


The choice of baseline L is particularly critical, in the sense that at the time that 
a neutrino factory would start running, there will be already existing experiments 
at baselines from 730 km from FNAL and CERN: (1) at Soudan, MINOS with its 
5.4 kton fiducial mass will have been fully operational, and similarly, (2) at GranSasso, 
the OPERA experiment and a multi-kton ICARUS detector. In view of the existence 
of such massive devices, it is worth considering if these detectors at their current 
baselines could be reused or improved in the context of the neutrino factory. If new 
sites have to be found in order to satisfy the requirements of longer baselines, major 
new “investments” will be required. 


As far as the neutrino energy is concerned, it is clear that the average neutrino 
energy E v scales linearly with the muon beam energy E^. A non-negligible aspect of 
the neutrino factory is the need to accelerate quickly the muons to the desired energy, 
and so, it is expected that higher energies will be more demanding that lower ones. 
Eventually, cost arguments could determine the muon energy. It could therefore be 
that lower energy, more intense neutrino factories could be more advantageous than 
higher energy, less intense ones. 


2 The effects related to 5 in vacuum 


In a three-family neutrino oscillation scenario, the flavor eigenstates v a (o. = e, /i, r) are 
related to the mass eigenstates v[{i = 1, 2, 3) by the mixing matrix U 


v n = U m i/ 


( 1 ) 


and it is customary to parameterize it as: 


C12C13 

U(9i2,9i3,9 2 3 , 5) = | —S12C23 — c 12Sl3 s 23 e 

S12S23 - c 12 s 13 c 2 3 e l5 


iS 


iS 


S12C13 

C12C23 — Si2Si 3 S 2 3e 

— C12S23 — Si 2 Si 3 C23e 


-iS 


iS 


S I 3 e 

Cl3 s 23 

Cl3 c 23 


( 2 ) 
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with Sij = sin % and c tJ = cos 0 l3 . 

We adopt the usual neutrino mass assignment in which the smallest = 

m\ — m{ is assigned to the solar neutrino deficit and the largest mass difference, 
A mh ~ Am | 2 is describing the atmospheric neutrino observations. We neglect ac¬ 
cordingly the important results from LSND, awaiting further clarification from the 
MiniBOONE experiment. We therefore consider the following indicative set of oscil¬ 
lation parameters^, compatible with current experimental observations (neglecting the 
LSND claim): 


Am 32 = 3 x 10 3 eV 2 , sin 2 d 2 3 = 0.5 (3) 

The value of the angle $13 is currently not known, and the best experimental bound 
comes from CHOOZ|T2[ which limits sin 2 26h 3 < 0.1. We take: 


sin 2 20 13 = 0.05 (4) 

For the solar parameters, we assume values compatible with the LMA-solution, 


which is known to yield optimal conditions for CP-violation studies (see Ref. [|T3| ): 


A-mli = 1 x 10 4 eV z , 


sin 2 di 2 = 0.5 


( 5 ) 


With this mass assignment and with the parameterization of the mixing matrix 
described above, a neutrino factory can provide precise information on the largest 
squared mass difference Am 2 3 ~ Am | 3 and on both angles 6*23 and C .3 • 

In the case where neutrinos (or antineutrinos) propagate in vacuum, the probability 
for flavor transition has a simple behavior. It can be written as 

P{v a ~>Uf,-E v ,L) = J2 (6) 

jk 

where E u is the neutrino energy and L the baseline between the source and the detector. 
The Jarlskog factor is J a pjk = UpjU'p k U* K -U a y.. For antineutrinos, we must replace 
U > U , i.e. Jaf 3 jk y Ja 0 jk Jaf 3 kj J[ 3 ajk- 

Since in general we have J a / 3 jk 7 ^ Japjki the probability for neutrinos P(u a —> 
E v , L) can be different from that for anti-neutrinos P{p a —>■ up]E u ,L), leading to 
the CP-violating effects. Since J a pjk = Jp a jki there are also T-violating effects. From 
the definition of J also follows that (CPT-invariance) 

P(z/ Q -4 up] E v , L) = P(Dp w u a ] E v , L) (7) 

1 Note that in studies of effects related to S , the signs of the Ato 2 ’s are important. 
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The oscillation probability can be made more explicit by separating real and imag¬ 
inary contributions of J a / 3 jk■ 

P(z,„ ~ S„, L) = £ Jatjj + Y. (j«we~‘ Am ‘ i,,2E + (8) 

3 j>k 

= 5 a p - 4 ^2 JafSjk sin 2 {A jk ) (9) 

j>k 

+4 ^2 ^mJaPik sin (A jk ) cos (A jk ) 

j>k 

Here 

Ajk = A m 2 jk ^r (10) 

and A 12 + A 2 3 + A 3i = 0. For anti-neutrinos, we must replace S mJ a pj k — > — S mJ a fSjk■ 
Hence, the first two terms of P form a CP-even term, and the last one is a CP-odd 
term. 

In the rest of the paper, we concentrate on transitions between electron and muon 
neutrinos. One can explicitly calculate the probabilities, inserting the elements of the 
chosen mixing matrix parameterization: 

P{v e -> Vp) = P(P li -*• v e ) = 

4c? 3 sin 2 A 23 s 2 2 s 2 3 s 23 + c i 2 ( sin2 Ai 3 s 2 3 s 23 + sin 2 A l2 s\ 2 (l - (l + s 2 3 ) s 23 )) 

—-c 2 3 sin(2#i 2 )si 3 sin(2# 23 ) cos 5 [cos 2A 13 — cos 2A 23 — 2 cos(2d 12 ) sin 2 A 12 ] 

X 

+ 2 C i 3 sin § sin(2fl 12 )s 13 sin(20 23 ) [sin2A 12 - sin2Ai 3 + sin2A 23 ] (11) 

This expression has been split in a first part independent from the phase 6, and in 
the two parts proportional respectively to cos 5 and sin 6. To obtain the probabilities 
for Up —> u e and P e —> P M , we must replace 5 —> —5, with the effect of changing 
sin 6 —» —sin 5 and cos S —> cos S. The term proportional to sin<5 is the CP- or 
T-violating term, while the cos S term equally modifies the probability for both CP- 
conjugate states. 

The behavior in the range where the energy is too large or the baseline too short 
to be sensitive to the smallest mass difference, i.e. |Am 2 2 | <C E v /L or |A 12 | <C 1, is 
found by A 32 —» 0 and A 33 —» A 23 . Both terms depending on cos <5 and sin S vanish, 
and we recover the well known probability 

P(z4 -> cO ^ sin 2 (26>i 3 )s 23 sin 2 (A 23 ) (12) 

Going back to the general case, the v e —> oscillation probability can be simplified 

by introducing values for 6\ 2 and 0 23 consistent with current experimental knowledge. 
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Indeed, for = 6*23 = 7 t/ 4 , we find 


P( u e u ii) — 2 Cl 3 \ c i 3 s ^ n ^12 + s i 3 2 sin A 12 sin(A 13 + A 2 3 — 5) + (13) 

sin h sin 2 A 12 + 2si 3 (sin 2 A 13 + sin 2 A 23 ) 


From this dependence, we see that a precise measurement of the v e —> oscillation 
probability can yield information of the 5-phase provided that the other oscillation 
parameters in the expression are known sufficiently accurately. We also note that the 
phase h changes the neutrino energy position of the maximum of the oscillation, because 
of the term sin(Ai 3 + A 23 — 5). 

Clearly, the dependence of the parameter h is a priori most “visible” in the energy- 
baseline range such that |A 12 | = | AE V ~ 1 and |A 23 | = \Am, 2 3 \L/AE u ~ 1. 
This is the region where all terms in Eq. (p~3|) contribute. 

When | A 12 1 -C 1 and |A 23 | — 1, we obtain 


P(z/ e -*■ Vp) ~ ^c 2 3 |c 2 3 A 2 2 + 2s 2 3 (sin 2 A 13 + sin 2 A 23 ) 


(14) 


+2A 12 si3 


sin(Ai 3 + A 23 ) cos h + (1 - cos(Ai 3 + A 23 )) sin 5 


At even higher E v or smaller L , we further have, when both |A 32 1 1 and 

A 13 1, |A 23 | <C 1: 


Piye 


V,,) ^ 


-13 


r 2 A 2 
'13^12 


+ 2s 2 3 (A^g + A 23 ) + 2 A 12 si 3 (A 13 + A 23 ) cosh 


and the dependence on the phase is only through cosh. From this follows a degeneracy 
under the change of sign of h, as was pointed out in Ref. |L 6 |. So, in this L and E u 


range, a precise determination of the oscillation probability can no longer determine 
the sign of h. 


The behavior for the two baselines L = 730 km and 2900 km as a function of 
neutrino energy E v is explicitly shown in Figure |lj. The probabilities are computed for 
three values of the h-phase: h = 0 (line), h = + 7 t /2 (dashed), h = — 7 t /2 (dotted). The 
other oscillation parameters are those described in Eqs. |3] and [5|: Am 32 = 3 x 10 ~ 3 eV 2 , 
Am 21 = 1 x 10 -4 eV 2 , sin 2 0 23 = 0.5, sin 2 9 12 = 0.5, and sin 2 29i 3 = 0.05. 

A region corresponding to the oscillation of the “first maximum” is clearly visible 
on the curves. We define the energy of the “first maximum” as follows 


7T 


A 32 — — 


jjjmax — 




jjjmax 


2tt 

(GeV) ~ Ani 3 2 (eV 2 ) 


2 x 1.27\ 


L(km) 


(16) 
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which yields E™ ax ~ 2 GeV at 730 km, E™ ax ~ 8 GeV at 2900 km and E™ ax ~ 20 GeV 
at 7400 km for A m\ 2 = 3 x 10~ 3 eV 2 . This energy corresponds to the point of maximum 
oscillation induced by Am| 2 and coincides with the maximum when <5 = 0. It will be 
useful when we discuss the point of maximum sensitivity to the 5-phase. 

A similar relation can also be defined in terms of L/E u which is 

L \ max 2tt f L(km) \ max _ ( 7r 

E^J = Am 2 2 " \E iy (GeV)) ~ \2 x 1.27Am| 2 (eV 2 ) 

which is approximately equal ( L/E u ) max ~ 400 krn/GeV for A m\ 2 = 3 x 10 -3 eV 2 . 

The dependence at high energy is not easily understood from a plot of the prob¬ 
abilities as a function of energy. In addition, while the effect is expected to be most 
“visible” at low energies, experiments at low energy are a priori more difficult to per¬ 
form than at high energy, because of the rising cross-section, easiness of detection and 
neutrino flux considerations. In fact, what is more important is not the probability 
behavior, but rather how will the number of events, observed in a given experiment 
at a given baseline L within a certain neutrino energy range, depend on the phase <5, 
relative to the total number of observed events in absence of 5 phase effects. Hence, in 
the following, we consider for definiteness the neutrino flux of a neutrino factory and 
will in order to ease the comparison between various energies and baselines introduce 
a “rescaled” probability. 



3 The neutrino factory (NF) 


The decay of the muon is a very well known process. Within the Standard Model and 
to leading order, the neutrino energy distribution in the muon rest frame in /i~ decays 
is given by 


d' 2 N 2r 2 

^<x-[(3-2,) + (l— 2,)/>„c„s0] 
d 2 N,j 12x 2 r/ . . ^ 


(18) 

(19) 


where is the average polarization of the muon beam, x = 2E lJ /rn IJ1 and 9 is the 
angle between the momentum vector of the neutrino and the mean angle of the muon 
polarization. The beam polarization can be carefully measured via the spectrum of the 
electrons from muon decay; thus the spectra of the two components of the neutrino 
beam can be known with good accuracy, and so the ratio of the fluxes reaching the far 
detector. 


Assuming perfect focusing, an unpolarized muon beam and for very relativistic 
muons, the energy spectrum of neutrinos detected at small angles with respect to the 
muon flight direction as in the case of long baseline experiment, is given by 

4>„„{z) oc 2z 2 (3 - 2z); 0 Pe (z ) oc 12^ 2 (1 - z) (20) 
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where z = E v jE^. A key point is that angular opening of the neutrino beam will shrink 
due to the Lorentz boost as y~ 2 = (E^/m^ 2 and hence the total neutrino fluency at 
a far detector will increase as E 2 . This means that the flux can be expressed as 

$«/„ (-£#*> z ) = z ) = jjM z ) (21) 


The scaling of the neutrino event rate can be expressed as 


N v 



(E^ L, E u ) a u (E u )dE u 


E 2 f l 

— J 4 > u (z)a u (zEf 1 )dz 


( 22 ) 

(23) 


since dE u = E^dz. By approximating the neutrino cross-section as a linear function of 
the neutrino energy, i.e. a v (zEfj) ~ crj) x E u , where is a constant, we obtain 


E 3 r 1 

N„ '-z J <!>zZdz OC -JT- 


(24) 


The integral over 0 is independent of E M and we find the known rapid growth of neutrino 
events with E^, characteristic of the neutrino factory. 


4 The rescaled probabilities 


In order to compare effects at different energies and various baselines, we define a 
“rescaled probability” parameter that allows a direct comparison of effects. 


Since to a good approximation (when z < 0.6 for electron-neutrinos and always for 
muon-neutrinos) the neutrino energy distribution behaves like E 2 (see section [|) and 
in addition, the neutrino flux scales like L~ 2 due to the beam divergence, we define the 
“rescaled probability” parameter p(u a up; E v , L ) as 


p(u a -4 up; E v , L) = P(u a -w up; E v , L) X 


(25) 


We note that 


1. it approximately correctly “weighs” the probability by the neutrino energy spec¬ 
trum El of the neutrino factory spectrum; 

2. it can be directly compared at different baselines, since it contains the L~ 2 at¬ 
tenuation of the neutrino flux with distance L ; 

3. p tends to a constant for E v —> 00 , hence the high energy behavior can be easily 
studied. 
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We start by illustrating the behavior of the rescaled probability for v e —► oscil¬ 
lations (see Eq. [15]) in the case where both |A 12 | <C 1 and |A 13 |, |A 23 1 -C 1. We note 
that the expression depends on the second power of the A^’s. The resulting ( L/E u ) 2 
dependence is canceled in the rescaled probability: 

p ( u e zy t ) ~ 2 C 1 3 { C 1 3 A 12+ 2 si 3 (A? 3 + A 23 ) + 2Ai 2 si 3 (Ai3 + A 23 )cos5j x 

- ^ c i 3 | c i 3 (A^-i 2 ) 2 + 2s 33 ((A m 2 13 ) 2 + (Am^ 3 ) 2 ) 

+8Am 32 Si 3 (Am 33 + Am 23 ) cos5j (26) 

As expected, the result is independent of the energy E v and the baseline L , and the 
rescaled probability “constant” is modified with a cos <5 dependent term. The sign of 
this term depends as expected on the sign of the Arne’s. 

Figure |2] shows the rescaled probability for v e —»■ oscillations in vacuum for three 
baselines L = 730 km (line), 2900 km (dashed) and 7400 km (dotted) as a function of 
neutrino energy E u . In this case, we chose 5 = 0 in order to illustrate the effect of the 
rescaling. The interesting oscillations at low energy driven by the interference of the 
A 32 , A i3 and A 23 terms are largely “damped” by the rescaling, and the situation has 
dramatically changed compared to Figure [Lj. 

We also note that the asymptotic value reached at high energy is independent of 
the distance L, reflecting the well-known result that the oscillated number of events 
is independent of distance L at high energy, since the L 2 growing of the oscillation 
cancels the L~ 2 attenuation of the flux with distance. 

The dependence on the phase 5 is illustrated in Figure [| For 5 = ±vr/2, the cos 5 
term vanishes and the probability is reduced at high energy compared to the 5 = 0 
case. The asymptotic value reached in the high energy limit is obviously independent 
of the baseline and is different than that reached for 5 = 0. For 5 = ±7t/2, it converges 
to the same values when the condition for |A 12 | <C 1 and |A 13 |, |A 23 | 1 is true. This 

depends on the baseline L considered. For L = 730(2900) km, the probabilities for 
5 = ±7t/2 are indistinguishable to within 10% for E u > 20(80) GeV. 


5 Propagation in matter 

Since neutrino factories will be associated to long baseline, it is not possible to avoid 
including effects associated to the neutrino propagation through the Earth matter. The 
simplest way to take into account these effects is to maintain the formalism developed 
for propagation in vacuum and to replace the mixing angles and the neutrino mass 
differences by “effective” values. 

The general exact and precise treatment of the effects induced by the propagation 
through the Earth can be quite complicated, due to the a priori complexity of the 


matter profile in the Earth. We will hence assume that the effects due to matter 
can be parameterized by one parameter, i.e. the “average” density p, to reproduce 
approximately the effect of traversing the Earth. This assumption was recently revisited 
in Ref. |18|], where it was concluded that for a baseline of 3000 km or less, the effect 
of the matter profile is not important. For very long baselines like L fa 7400 km, 
they conclude that the Earth profile is important. This however does not pose us a 
problem, since, as we will show below, we will not consider that baseline in the context 
of the searches for 5-phase-induced effects. We hence all along use the approximation 
of constant density. 


In the case of two-neutrino oscillation, in which the mixing is driven by a single 
mixing angle 9, and propagating in matter of constant density p, it is rather easy to 
understand the behavior of the oscillation probability for neutrinos or antineutrinos. 
An important quantity for matter effects is D, defined as 


D(E v ,p) = 2V2G F n e E v = 7.56 x = -D(E Pl p) (27) 

gem 6 Gev 


where n e is the electron density and p the matter density. For antineutrinos, D is 
replaced by — D. 


I 11 the case of three neutrino mixing, the situation is mathematically more complex. 
Approximate oscillation probabilities have been considered by perturbative expansion 
in the oscillation parameters |T9|, [20| , [L5|]. Exact formulas have been derived^], |23j 
and we make use of this formalism as implemented in Ref. [17]]. For completeness, we 


reproduce here the expressions. The mass eigenvalues in matter M x , M 2 and M 3 are: 


Ml = m 2 + - - -VA 2 - 3 BS - —VA 2 - ZBVl^^S 2 
3 3 3 

Ml = m? + — — -VA 2 - 3 BS + —VA 2 - WV^S 2 
2 1 3 3 3 

Ml = m? H-1— \/A 2 — 3BS 

3 3 


(28) 

(29) 

(30) 


where A, B and S are given in the Appendix of Ref. 
matter the analytical expressions read: 


17| . For the mixing angles in 


sin 2 0™ 
sin 2 0™ 
sin 2 0 2 m s 


-(M 2 4 - aM'l + /3) A Ml 
AM| 2 (M 4 - aM 2 + (3) — AM| 1 (M| - aM 2 + (3) 
M 4 - aM 2 + (3 
AMIAM 2 , 

G ,2g 23 + F 2<: 23 + 2GFC23S23CS 

G 2 + F 2 


(31) 

(32) 

(33) 


where a, f3 , G and F are found in the Appendix of Ref. 
for the 5-phase in matter, i.e. S m . 


T7|, and also an expression 
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It can be shown that in the relevant situation in which (Am^l <C (Am^l ~ |Am| 2 |, 
matter effects in three-family mixing decouple to two independent two-family mix¬ 
ing scenarios. This is essentially because matter effects become important when the 
parameter D(E u ,p) is similar to one of the |Am 2 j's and hence the matter effects 
driven by (Am^l will occur for fixed density at a very different energy E u than for 
Am 2 l rs |At? 7 J 2 |. The term |A? 7 i 21 | (resp. jAm 2 ^) will produce a MSW resonance in 
the #12 (resp. # 13 ) angle. 


Given the energies and baselines considered, we illustrate analytically our arguments 
for a two-family mixing matter effect in v e —> v M oscillations, in which we identify the 
2-mixing angle # with the 3-mixing # 13 and the mass difference squared as (Am^l ss 
|Am 2 2 |. In the case |A 12 | <C |A 13 | ~ |A 


23 1 


#™, #™ and 6" 


1 in matter, the effective mixing angles 
will be approximately equal to (see e.g. Figure 1 of Ref. ||I7|P 2 


sin 2 6™(D) ~ 

sin 2 2 #™(£>) -> 

sin 2 2 #™(L>) = 


5 m « 5 


sin 2 #23 
Ri 0 

sin 2 2# 13 

sin 2 2#i 3 + (^|- - cos 2 #i 3 ) 2 

sin 2 2 #i 3 
sin 2 2# 13 + M 2 


(34) 


where A4 = D / Am\ 2 — 1 and, for the second line, we assume # 13 small. This term 
has the well-known MSW resonance behavior, which predicts that the effective angle 
sin 2 2 0™z(D) —> 1 for D —> Am 2 2 (i.e. A4 —> 0), no matter how small sin 2 2# 13 is. The 
effective mass squared difference becomes 

AM 32 ~ Am 2 2 \/sin 2 2# 13 + A4 2 (35) 


which implies that the effective oscillation wavelength becomes very large (resp. small) 
close to (resp. far from) the resonance energy. 

The oscillation probability at small L/E v (see Eq. [L2D is in matter 
P m (ve^ v») - sin 2 (2#-)4sin 2 (A- 3 ) 


sin 2 2 # 


13 


sin 2 2# 13 + M 


sin sin" 


2 15 23 


2 ( Am\ 2 \/ sin 2 2# 13 + M 2 




4 E„ 


(36) 


In the case where the neutrino energy is at the MSW resonance (i.e. JA — 0 and 
E v = S res cos 2# 13 Am| 2 ), the probability is approximately 

P m (is e —»■ v u ) ~ So 3 sin 2 ( A ml 2 y/ sin 2 2# 13 —- 

_ v e 23 V 4^ es cos2# 13 Am 2 2 y 

2 In reality, the angles sin 2 0^ and sin 2 9^ tend to rise slightly for D > Am^, because the non¬ 
vanishing A rri 2 i splitting removes the degeneracy between muon and tau flavors (see Ref. {T|), but 
we neglect this effect. 
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3 23 


sin 


/ tan 26>13 \ 

V A£ res ) 


2 ( tan2g 13 y 2 

^23 l /^gres J 


( 37 ) 


where in the last line we assumed L < tan2#i 3 ,/4£ res « 14000 km. In this case, the 
probability grows with L 2 , giving rise to the enhanced oscillation probability at large 
distances through matter. 

For neutrino energies E v above the resonance energy El es , where one can approxi¬ 
mate M 2 > 1 > sin 2 20 13 , the probability is 

P m (u e ^U /1 ) ~ Sm ; j 2 gl3 g 23 sil12 ( 38 ) 

and this expression will vanish for growing A4 since the sine function will not cancel 
the M 2 at the denominator. These behaviors will be rediscussed later in the context 
of the 5-phase dependent terms. 

To summarize, there will be two specific neutrino energies of interest when neutrinos 
propagate through matter: 


1. for D ~ Am§ 2 , we reach for neutrinos the MSW resonance, in which the effective 
mixing angle sin 2 (2#^) ~ 1. In terms of neutrino energy, this implies 


= cos2» 13 Aml 2 =£TOcos 


2 y/2G F n e 
1.32 x 10 4 cos20 13 Am| 2 (eV 2 
p(g/cm 3 ) 


32 


in GeV 


(39) 


where £ res = (2\/2 Gpn e ) 1 . For density parameters p equal to 2.7, 3.2 and 
3.7 g/cm 3 one hnds El es ~ 14.1, 12.3 and 10.7 GeV for Am 2 2 = 3 x 10 ~ 3 eV 2 . 

2. for D > 2Aml 2 , the effective mixing angle for neutrinos is always smaller than 
that in vacuum, i.e. sin 2 (26*^) < sin 2 (2d 13 ). In terms of neutrino energy, this is 
equivalent to 

E v > 2 E r u es (40) 


3. these arguments are independent of the baseline L and depend only on the matter 
density p. 


In the rest of the paper, we will solve the equations numerically in order to properly 
introduce exact matters effects into our oscillation probabilities. 

The distortion of the oscillation probabilities introduced by the propagation through 
matter are shown in Figures [| and ^ where the probability for v e —> oscillations 
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for two baselines L = 730 km and 2900 km are shown as a function of neutrino energy 
E v . For each baseline, the probabilities are computed for three values of the h-phase: 
5 = 0, 5 = +7t/2, 5 = —7t/2. 

We further illustrate the behavior in matter in Figures | 6 |. |7] and |], where the rescaled 
probability for v e — > z/ /t oscillations for three baselines L = 730 km, 2900 km and 
7400 km are shown as a function of neutrino energy E v . For each baseline, the proba¬ 
bilities are computed for neutrinos in matter, in vacuum (dotted line, same for neutrinos 
and antineutrinos) and for antineutrinos in matter. A phase <5 = 0 has been assumed. 

The characteristic behavior in matter can be readily seen. At the shortest baseline, 
L = 730 km, the effect of matter is small, i.e. the probabilities are very similar to that 
in vacuum and the probability is simply slightly enhanced (suppressed) for neutrinos 
(antineutrinos). Between approximately 1 < E u < 10 GeV, the rescaled probability 
rises; this corresponds to the oscillation at the “first maximum” (see equation pT|) , 
which yields E™ ax ~ 2 GeV for 3 x 10 ~ 3 eV 2 . Above E v ~ 10 GeV, the rescaled 
probability reaches a constant. 

At the middle baseline L = 2900 km, the effect of propagation through matter 
is quite visible. The resonance for neutrinos can be seen at E xes (GeV ) ~ 12 GeV, 
which is quite close from the “first maximum”, which is now given the baseline at 
E™ ax ~ 3 GeV (see equation 0 )- As expected, the probability in vacuum reaches at 
higher energies the same constant as for the shortest baseline L = 730 km. However, 
in matter for E u > 2E res , the probability becomes smaller than that in vacuum. For 
antineutrinos, the probability tends to be suppressed by matter, but tends to the same 
value as that of neutrinos at high energy. 

For the longest baseline L = 7400 km, the effect of matter is quite strong. The “first 
maximum” is now at E™ ax ~ 20 GeV (see equation 0), but this corresponds to an 
energy very close to the condition for D > 2Am 2 2 , in which the effective mixing angle 
for neutrinos is smaller than that in vacuum. The probability at such large distances 
is therefore quite “suppressed” by matter as opposed to what it would be in vacuum. 

From these discussions, it is already clear that if one wants to study oscillations 
in the region of the “first maximum”, one should not choose a too large baseline L, 
otherwise, matter effects will suppress the oscillation probability. This will be further 
discussed in the next section. 


6 Detecting the S phase at the NF and the L/E u 
scaling 

We are now ready to investigate methods to optimally look for effects related to the 
phase 5 in neutrino oscillations. We consider that propagation of neutrinos will always 
occur through matter and hence use the exact numerical calculations for three family 
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mixing, without any approximation in the oscillation probabilities, nor in the treatment 
of matter effects. 

In order to further study the dependence of the 5-phase, we consider the following 
three quantities which are good discriminators for a non-vanishing phase 5: 

1. A 5 = P(u e —>• 5 = +7 t/2 ) - P{v e -> 1^,6 = 0) 

The discriminant A <5 can be used in an experiment where one is comparing the 
measured u e —> u fx oscillation probability as a function of the neutrino energy E u 
compared to a “Monte-Carlo prediction” of the spectrum in absence of 5-phase. 

2 . A C p(5) = P(u e -> v hi 5) - P(h e -> P M , 5) 

The discriminant A cp can be used in an experiment by comparing the appearance 
of // /t (resp. v M ) in a beam of stored /i + (resp. /j,~) decays as a function of the 
neutrino energy E v . 

3. A T (5) = P(y e -> v hi 8)-P(uh -> r'e, 5) or A T (5) = P(h e -> 5) -> P e , 5) 

The discriminant Ap can be used in an experiment by comparing the appearance 
of i/p (resp. P /t ) and P e (resp. u e ) and in a beam of stored /i + (resp. n~) decays 
as a function of the neutrino energy E v . 

Each of these discriminants have their advantages and disadvantages. 

The A5-method consists in searching for distortions in the visible energy spectrum 
of events produced by the 5-phase. While this method can in principle provide excellent 
determination of the phase limited only by the statistics of accumulated events, in 
practice, systematic effects will have to be carefully kept under control in order to look 
for a small effect in a seen-data versus Monte-Carlo-expected comparison. In addition, 
the precise knowledge of the other oscillation parameters will be important, and as 
will be discussed below, there is a risk of degeneracy between solutions and a possible 
strong correlation with the # 13 angle at high energy. 

The A C p is quite straight-forward, since it involves comparing the appearance of 
so-called wrong-sign muons for two polarities of the stored muon beam. It really takes 
advantage from the fact that experimentally energetic muons are rather easy to detect 
and identify due to their penetrating nature, and with the help of a magnetic field, 
their charge can be easily measured, in order to suppress the non-oscillated background 
from the beam. A non-vanishing A^p should in principle be a direct proof for a non¬ 
vanishing 5-phase. This method suffers, however, from the inability to perform long- 
baseline experiment through vacuum. Indeed, matter effects will largely “spoil” A cp 
since it involves both neutrinos and antineutrinos, which will oscillate very differently 
through matter. Hence, the A cp requires a good understanding of the effects related 
to matter. In addition, it involves measuring neutrinos and antineutrinos. The matter 
suppression of the antineutrinos will in practice determine the statistical accuracy with 
which the discriminant can be measured. 
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Finally, the Ap is the theoretically cleanest method, since it does not suffer from the 
problems of Ad and Aqp- Indeed, a difference in oscillation probabilities between v e —> 
vp and —»• v e would be a direct proof for a non-vanishing 5-phase. In addition, matter 
affects both probabilities in a same way, since it involves only neutrinos. Unfortunately, 
it is experimentally very challenging to discriminate the electron charge produced in 
the events, needed in order to suppress the background from the beam. However, one 
can decide to measure only neutrinos, which are enhanced by matter effects, as opposed 
to antineutrinos in the A cp which were matter suppressed, and hence the statistical 
accuracy of the measurement will be determined by the efficiency to recognize the 
electron charge, rather than by matter suppression. 


6.1 The L/E„ scaling 


Regardless of their advantages and disadvantages, there is one thing in common be¬ 
tween the three discriminants A 5 , A cp and A p: their behavior with respect to the 
neutrino energy E u and the baseline L. By explicit calculation, we find 


A 5 


--c ^ 3 sin 2 di 2 Si 3 sin 2023 X 


(41) 


j^COS 2 A13 
1 2 

_ 2 Cl 3 S 13 


— cos 2 A 2 3 — 2 cos 20 i 2 sin 2 Ai 2 + sin 2 Ai 2 — sin 2 A 13 + sin 2 A 2 3 
cos 2A 13 — cos 2A 23 + sin 2A 12 — sin 2 A 13 + sin 2A 23 ] 


where for the second line we assumed for simplicity 0i 2 = 0 2 3 = 7 t/ 4 , and similarly, 


A cp — A t 


c 2 3 s 13 sin 20 i 2 sin 20 23 
—c 2 3 Si 3 sin 2012 sin 20 


sin 5 


sin 2 A 


12 


23 


sin 5 


sin A 


12 


— sin 2 Ai 3 + sin 2 A 
sin A 13 sin A 23 ] 


23 


(42) 


As expected, both expressions vanish in the limit Am 2 2 —> 0 where Amf 3 —> Am | 3 
Also, as one reaches the higher energies, the terms A cp = A p vanish as 


|A( 7 p| = |Ap| ~ c 2 3 Si 3 sin 2012 sin 2023 sin 5Am 2 2 


4 E v 


sin 2 ( Amn 


23 


4 E v 


— c i 3 s i 3 s i n 2012 sin 20 2 3 sin 5Am\ 2 (Am 


2 \2 
23 / 


4A„ 


(43) 


hence, in the very high energy limit at fixed baseline, the effects decrease as E~ 3 . That 
the effects disappear at high energy is expected, since in this regime, the “oscillations” 
of the various A jfc ’s wash out. 

At high energy, we can also approximate cos2Ai 3 « cos 2 A 2 3 ~ 1, and find 

|A*| ~ ^A C p{d = 7t/ 2) = ^A t (5 = 7 t/2) (44) 
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so the CP- or T-conjugation is equivalent to a change of phase 5 —> —5, i.e. 5 = 
+7t/2 — ■> 5 = — 7t/2. 

The important point is that all expressions depend upon some factors which contain 
the various mixing angles, multiplied by oscillatory terms which always vary like sine 
or cosine of A jfc -terms (the terms in squared brackets in the expressions above). 

Hence, we expect the various discriminants to scale like Ajj. cx L/E v . The sensi¬ 
tivity to the 5-phase will therefore follow the behavior of the oscillation probability, 
and we therefore argue that the maximum of the effect will occur around the “first 
maximum” of the oscillations, i.e. for E™ ax = Am| 2 T/27r (see Eq. (|I6|)). 

Strictly speaking, the maximum of the 5-phase sensitivity does not lie exactly at 
the “first maximum” as defined in Eq. (|T6|). From Eq. (|4^), we expect the maximum 
to be “shifted” to higher values of L/E v , since it corresponds to the maximum of the 
term 

sin A 12 sin A i3 sin A 23 — A m i2^r sin2 (^ m 23 (45) 

which has the functional form x sin 2 x and, therefore, has its maximum shifted to 
higher values of x compared to sin 2 x. This small shift is smaller than the oscillation 
wavelength itself, and does not cause a major problem, since experimentally we will 
always have sufficient energy range to cover the full oscillation. 

These considerations are strictly true only for propagation in vacuum. When neu¬ 
trinos propagate through matter, matter effects will alter these conclusions. We will 
however show that as long as the baseline is smaller than some distance such that the 
corresponding “first maximum” E™ ax lies below the MSW resonance neutrino energy 
E xes , the considerations related to the L/E u scaling are still largely valid. 

To illustrate this, we show in Figure || the rescaled probability for v e —> v M oscilla¬ 
tions involving neutrinos and as a function of L/E v , for three baselines L = 730 km, 
2900 km, 7400 km and in vacuum (independent of baseline). We assume 5 = 0. We 
readily observe that for the shortest baseline L = 730 km the curve is very similar to 
the one in vacuum. For the middle baseline L = 2900 km, the curve is modified with 
respect to the vacuum case, in the sense that the probability around ( L/E u ) max ~ 400 
(see Eq.[T7]) is enhanced, but for ( L/E v ) < 130 the probability is smaller than that in 
vacuum. For the longest baseline L = 7400 km, the curve is highly distorted by matter 
and already for (L/E u ) < 400 is the probability smaller than that in vacuum. 

In what way does the matter effect alter the ability to look for effects related to 
the 5-phase? It is incorrect to believe that only the A cp discriminant will be affected 
by propagation through matter, since it is the only one to a priori mix neutrinos and 
antineutrinos. In reality, the “dangerous” effect of matter is to reduce the dependence 
of the probability on the 5-phase, and this for any kind of discriminant. 

I 11 matter, we would for example write the A t discriminant as 


|A™I 




3 m 
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: 13 Rs™ sin 20(^ sin 20 2 3 sin A 


(c™i 2 - m 


sin A™ sin A™ sin A™ 


(46) 


where because of our choice of Am^’s 


we have 6 ^ 


0 2 3 and A m « A (see Eq. (|34|)). 
This implies that the A-phase discriminants have a different structure that the terms 
that define the probability of the oscillation (i.e. the non A-phase dependent terms). 
The discriminants are the products of sines and cosines of all mixing angles and of the 
Ajfc’s (see Eqs. (flTD and 
different. 


|)). Because of this structure, their property in matter is 


The Jarlskog’s determinant can be shown to be “invariant” under matter effects|25 
in the sense that the product 


J ^Am| 2 Am 2 3 Am 2 3 


= CJ 3 S 13 sin 2 0 12 sin 20 2 3 sin A ( Amj 2 Amj 3 Am; 


"23 


(47) 


is independent from matter effects. This point was also discussed in Ref. |T(| , where it 
was pointed out that this implies that the A-phase terms will be independent from the 
matter effects as long as we can approximate 


sin A 12 sin A 13 sin A 


*23 


A12A13A23. 


(48) 


This fact is not very revelant in the current context, since the best sensitivity to the 
A-phase is expected just when this approximation is not valid (i.e. |A 13 | ~ | A 2 3 1 ~ 1), 
otherwise the effects wash out. As a consequence, we expect that the CP-odd terms 
do depend on matter effects in the relevant baseline and energy region. 

The behavior for neutrino energies above the MSW resonance ff” es is determined 
by the fact that in this energy regime, 0^(E > E” es ) —> 7 t/ 2 and therefore —> 0 and 
s i 3 ~* 1 - More explicitly, one can show that 


m\2 m 
l 13/ *13 

oc 

hi 
* 1 

to 

sin 2 6^ 


sin 2023 = const. 

sin 2 0™ 

-► 

const. 

CNCO 
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A M| 2 oc E v 

A Ml 


A m 32 = const 


Therefore, 


|A!p| 


oc E„ 2 sin A 


sin ( A 


32 


w) sil1 ' 2 ( AM « 


4 E v 


oc E, 


-3 


(49) 


(50) 


and one recovers a neutrino energy dependence identical to that in vacuum (see Eq. (|43|) ) 
Note also that the argument of the sine function AM 2 3 L/4E„ is not small (i.e. the ap¬ 
proximation sin x ~ x is not valid). For our choice of oscillation parameters, the mass 
difference is approximately equal to AM 2 3 (eV 2 ) ~ 3 x 10 ” 4 x E I/ (GeV), and hence the 
dependence on the baseline is 


| Ay | 


oc E„ 2 sin A 


sm 




Sill" 


(3.8 x 10” 4 L(km)) (51) 
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Hence, the discriminant will first be enhanced and then be suppressed by matter effects. 
The maximum is found when the sine squared function reaches a maximum, or at 
approximately 4000 km under the assumption of high energy neutrinos. We will verify 
numerically that this is really the case in Section ^]. 

As anticipated, these discussions say that if one wants to study oscillations in the 
region of the “first maximum”, one should not choose a too large baseline L, otherwise, 
matter effects will suppress the oscillation probability. This is even more so true, as 
it will be recalled below, that the magnitude of the effects related to the 5-phase are 
suppressed more rapidly than the oscillation. 

The simplest way to express the condition on the matter is to require that the 
energy of the “first maximum” be smaller than the MSW resonance energy: 

2V2G F n e E™ ax < Am\ 2 cos 2 0 13 (52) 


and, by inserting the 
L < 

J - / max rv 


definition of E™ ax = Am 32 L/2ir we get 
7 r cos 26*13 7 r cos 26*13 


V2G F n e 2 x 1.27 x 7.56 x 10 - 5 (eV 2 )p(g/cm 3 ) 
1.5 x 10 4 (km) 


p(g/cm 3 ) 


5000 km 


(53) 


To summarize, we find that the discriminants of the 6-phase all scale with L/E v . 
The most favorable choice of neutrino energy E„ and baseline L is in the region of 
the “first maximum” given by ( L/E u ) max ~ 400 for |Am§ 2 | = 3 x 10 ~ 3 eV 2 . This 
leaves a great flexibility in the choice of the actual neutrino energy and the baseline, 
since only their ratio L/E u is determinant. Because of the rising neutrino cross-section 
with energy, we will see below that it will more favorable to go to higher energies 
if the neutrino fluency is constant. Keeping the L/E v ratio constant, this implies an 
optimization at longer baselines L. One will hence gain with the baseline L until we 
reach L max pc 5000 km beyond which matter effects will spoil our sensitivity. 


7 Detection of at the NF 

We begin the discussion with the detection of the 5-phase with the method of the 
A 5 discriminant. We recall that this method implies the comparison of the measured 
iz e —> oscillation probability as a function of the neutrino energy E v compared to a 
“Monte-Carlo prediction” of the spectrum in absence of 5-phase. 

We first illustrate that the discriminant really scales like L/E v in Figure 0 . where 
the rescaled A $ discriminant is shown as a function of the L/E u ratio. The plot is 
computed for neutrinos propagating in matter at three different baselines L = 730 km, 
2900 km and 7400 km, and for neutrinos propagating in vacuum. Indeed, for the 
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shortest two baselines L = 730 km and 2900 km, the value of the discriminant is 
very similar to that obtained in vacuum. For the longest baseline L = 7400 km, well 
above L max (see Eq. (|53|) ), the effect is highly suppressed by matter effect, since the 
“first maximum” corresponds for the given baseline L to an energy E u above the MSW 
resonance energy E r u es . 

From the figure, we clearly also see that the effect is largest around the “first 
maximum”, in the region given by ( L/E v ) max ~ 400. 

Given the freedom in the energy E v and baseline L, provided the ratio L/E u is 
appropriately chosen, and given the matter effects at distances beyond ~ 5000 km, we 
restrict our choices to the two baselines L = 730 km and 2900 km. 


7.1 The correlation with #13 


When searching for effects related to the 5-phase by comparing the measured v e —> 
oscillation probability as a function of the neutrino energy E v to a “Monte-Carlo pre¬ 
diction” of the spectrum in absence of 5-phase, requires necessarily a precise knowledge 
of the other oscillation parameters entering in the oscillation probability expression. 

In particular, the knowledge of the angle #13 could be quite important. Indeed, the 
v e —> v M oscillation is primarily driven by the 9 13 angle and only thanks to a different 
energy dependence of the terms proportional to 5 than to those independent of 5 can 
one hope to determine 9 13 and 5 at the same time! 

This is however not true at high energy, when both IA 12 I -C 1 and |Ai 3 |, |A 23 1 -C 1. 
This can be explicitly shown for example for simplicity in the limit of small 0 13 . The 
rescaled probability is in this case (see Eq. (|26|) ) a constant: 

piy e -> cO — — 7^— { 1 + 2M4 + 8Ns 13 cos 51 (54) 

where M = ((Am ^) 2 + (Am| 3 ) 2 )/(Amf 2 ) 2 and N = (Am 2 3 + A/n| 3 )/(Am 2 2 ). The 
absence of “oscillations” at high energy implies that a change of d 13 can mimic a change 
of 5. 


To illustrate this effect in a concrete example, we show in Figures 0 and [ 12 ] the 
rescaled probabilities for the two baselines L = 730 km and 2900 km as a function of 
neutrino energy E v . The probabilities are computed for neutrinos in matter (full line) 
and in vacuum (dotted line), and for three values of the 5-phase: 5 = 0, 5 = +7 t/2 , 
5 = — 7t/2. 


At the shortest baseline L = 730 km, this behavior can be clearly seen for E u > 
20 GeV. The two curves for 5 = + 7 t /2 and 5 = —tt/2 tend to a constant, which is 
clearly different than the constant for 5 = 0, but that cannot be distinguished from a 
5 = 0 with a different 9 13. In addition, we recall that the reason that both 5 = + 7 t /2 
and 5 = —tt/2 tend to the same constant, is due to the fact that the probability at 
high energy depends only on cos 5 (see Eq. (|26|)). 


18 




At the longer baseline L = 2900 km, the situation is improved since up to ~ 
100 GeV, the rescaled probability is clearly not a constant, but rather still falling with 
energy. We also point out that one can visually see that most of the sensitivity at 
L = 2900 km is for neutrino energies around 10 GeV or so. This is an important 
point to take into account when considering the muon detection threshold of a given 
experimental setup. 


8 Direct detection of T -violation at the NF 

We now turn to the most challenging search for effects related to 5 and consider the A t 
discriminant. We recall here that the discriminant At implies the comparison between 
the appearance of (resp. and v e (resp. u e ) in a beam of stored /r + (resp. n~) 
decays as a function of the neutrino energy E v . 

The A t method takes in advantage that deals with differences between neutrinos 
(or anti-neutrinos) only, so matter effects give only a scale factor and no “CP-fake” 
violation effects are expected for 5=0. So, the comparison of v,, —> u e and u e —> 
oscillation probabilities offers a direct way to highlight a complex component in the 
mixing matrix, independent of matter and other oscillation parameters. 

This is clearly visible in Figure [13], where the rescaled value of A t is shown for 
three baselines and for vacuum. The magnitude of the T-violation effect is roughly 
the same for vacuum, 730 km and 2900 km, while for very long baselines (7400 km) 
the sensitivity is lost. As described in Section |6.1|, all 5-phase discriminants scale with 
L/E u until 1^5000 km (see equation [53]) and because of matter effects, beyond this 
value the effect decreases dramatically. 

The drop of the rescaled A r for very long baselines is also illustrated in Figure [14]. 
The rescaled At discriminant is shown as a function of the baseline and of the neutrino 
energy, for neutrinos in vacuum (top) and in matter (bottom) for 5 = +7 t/2 . In vacuum, 
the same value of At x E 2 / L 2 is obtained keeping the L/E v ratio constant, while in 
matter this pattern is destroyed, loosing in sensitivity for long baselines. For instance, 
for E v ~ 12 GeV, the discriminant starts dropping for baselines larger than 4000 km. 

As already pointed out, the major difficulty of the A T method is that it requires the 
measurement of the electron charge in order to discriminate the large v e background 
from the beam. Both, the wrong-sign muon and the wrong-sign electron samples are 
needed to build the magnitude. 


9 Direct detection of CP-violation at the NF 

We now continue our discussion and consider the A cp discriminant. We recall that the 
A cp is in principle quite straight-forward, since it involves comparing the appearance 
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of Up (resp. Pfj) in a beam of stored /i + (resp. p - ) decays as a function of the neutrino 
energy E v . A non-vanishing A cp obtained from an observation of the two polarities of 
the stored muon beam would be a direct proof for a non-vanishing 5-phase. 

However, as already pointed out, this method suffers from the practical necessity 
to traverse the Earth matter in order to perform long baseline experiments. 

Since neutrinos and anti-neutrinos are affected by matter with an “opposite sign”, 
the matter terms will not cancel in the A cp discriminant. As a result, the A cp will 
not necessarily vanish when 5 = 0 (having the so called “CP-fake” violation effects. 
Hence, a direct detection of a non-vanishing A cp discriminant does not imply 5^0. 

We illustrate this situation in Figure [I5|, where the rescaled A cp discriminant is 
plotted as a function of the neutrino energy E u , computed for neutrinos propagating 
in matter at three different baselines L = 730 km, 2900 km and 7400 km. Three sets 
of curves are represented, corresponding to 5 = + 7 t/ 2 , 5 = —tt/2 and 5 = 0. 

At the shortest baseline L = 730 km, we observe that the maximum of the effect 
is indeed around the “first maximum” at about a neutrino energy of 2 GeV. For 5 = 0 
(thin curve), the discriminant does not vanish, due to matter. This effect is however 
smaller than the “genuine” effect at |5| = 7 t /2 and we observe that for 5 = + 7 t /2 —> 
— 7 t /2 the discriminant changes sign. 

For the medium baseline of L — 2900 km, the maximum of the effect is again 
near the “first maximum” (E u = 8 GeV). The A cp x E 2 / L 2 discriminant is larger (a 
factor 2) than for L = 730 km in the case of 5 = + 7 t/ 2 , but this includes the “CP- 
fake” contribution coming from matter. In order to compare at both baselines the 
“genuine” CP-violation effect (the one due to a non zero 5-phase) one has to compute 
the difference between 5 = + 7 t /2 and 5 = —vr/2, resulting similar at both baselines. 
Another interesting feature is that, at L — 2900 km the discriminant is already positive 
for any value of 5. 

At larger baselines (L = 7400 km) A C p is completely dominated by matter effects. 
The three values of 5 mix and the sensitivity on A^p to the “genuine” CP -violation 
is lost. 

Finally, we address the scaling of the effect with LIE u . Figure [16] shows the same set 
of curve displayed in Figure [L5|, but as a function of the L/E v ratio. In addition to the 
three different baselines, the case of neutrinos propagating in vacuum (independent of 
baseline) is also shown for comparison. At 730 km and 2900 km the effect is maximum 
for the same L/E v ~ 400 km/GeV (for A= 3 x 10 ~ 3 eV 2 ). In addition, the 
difference between 5 = + 7 t /2 and — 7 t /2 at the maximum is similar for both baselines. 
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10 Behavior for small Q\% values 


In the above sections, the possibility of detecting CP-violation, based on the compar¬ 
ison of oscillations involving electron and muon (anti-)neutrinos has been discussed. 
This transition is strongly dependent on the value of the mixing angle # 13 , and it is 
natural to ask how the sensitivity to CP-violation varies with the value of this angle. 
It was stated (|26|]) that, “as far as the Am 21 effects in the CP-conserving part of 
the oscillation probability can be neglected”, the sensitivity to CP-violation does not 
depend on 0 13 . 

We can check this hypothesis directly with the help of our Eq. 0 . The probability 
for oscillations depends on 


4c 2 


sin 2 A 23 Si2 s i3 s 23 + C 12 ( sin2 Ai3S? 3 s 23 + sin 2 A i 2 s? 2 ( 1 - ( 1 + s? 3 ) s\ 


’23 


sin A 23 s 12 £l 3 S 23 + c 12 (sin- 
where we have neglected A 12 . The parts dependent on 5 go like 

— -c 2 3 sin( 2 # 12 )s 13 sin(20 23 ) cos 5 [cos2Ai 3 — cos 2 A 23 — 2 cos(2$i 2 ) sin 2 A 12 ] 

+ ^c 2 3 sin 5 sin(26> i2 )gi 3 sin(26> 23 ) [sin 2 Ai 2 - sin 2 Ai 3 + sin 2 A 23 ] 


(55) 


(56) 


As can be seen, the 5-dependent part of the probability depends on sin C 3 , while 
the leading term (and the total number of oscillated events) is proportional to sin 2 # 13 . 

In the Gaussian approximation, the statistical error is proportional to the square 
root of the number of events, so the significance of the 5-effects is independent of C 3 . 
The validity of this assumption will depend on the total number of muon decays and 
on the total mass of the considered detector. 


As far as the Am 21 part is concerned, the approximation is valid if 

sin 2 Aiss^s^ > sin 2 A 12 s 2 2 (l - (l + £ 13 ) s 23 ) (57) 
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(58) 


and Am | 3 = 3 x 10~ 3 eV 2 . For smaller values of d 13 , the oscillation probability is 
dominated by the constant term, and the cancellation of the # 13 dependence no longer 
applies. 


11 Two concrete cases at 730 and 2900 km 

In order to assess with concrete examples the use of the 5-phase discriminants, we 
consider the two baselines, with corresponding muon beam energy and matter densities: 
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Baseline 

E, 

Muon decays 

Matter density 

732 km 

2900 km 

7.5 GeV 
30 GeV 

10 21 

2.5 x 10 20 

2.8 g/cm 3 

3.2 g/cm 3 


Both examples were chosen to have the same L/E M and include matter effects. 
Because of the linear rise of the neutrino cross-section with E u , the factor 4 in muon 
energy between the 732 km and 2900 km case, is “compensated” by an increase of 
intensity by the same factor in favor of the shorter baseline. 

We compute the fluxes assuming unpolarized muons and disregarding muon beams 
divergences within the storage ring. We consider that 50% of the stored muons decay 
on the direction of the detector (“useful” muons). We integrate the expected event 
rates using a neutrino-nucleon Monte-Carlo generator [^7| . The total charged current 
(CC) cross section is technically subdivided into three parts: the exclusive quasi-elastic 
scattering channel oqe and the inelastic cross section a ine i asic which includes all other 
processes except charm production which is included separately. Event rates for the 
two energy-baselinc-flux combinations are shown in Tables [L] and |2| for a 10 kton fiducial 
mass detector. 


Our analyses are performed on samples of fully generated Monte-Carlo events p7| 


which include proper kinematics of the events, full hadronization of the recoiling jet 
and proper exclusive polarized tau decays when relevant. Nuclear effects, which are 
taken into account by the FLUKA model, are included as they are important for a 
proper estimation of the tau kinematical identification. 


To have a proper account for the detector and background effects, the specific 
detector model of a Liquid Argon TPC|| has been assumed. 

The detector response is included in our analyses using a fast simulation which 
parameterizes the momentum and angular resolution of the emerging particles, using 
essentially the following values: electromagnetic shower 3%/^/EQl%, hadronic shower 
~ 20 %/\fE © 5%, and magnetic muon momentum measurement 20%. 


11.1 Lepton charge identification 

Muon identification, charge and momentum measurement provide discrimination be¬ 
tween and charged current (CC) events. Good v e CC versus v NC discrimination 
relies on the fine granularity of the target. Finally, the identification of v T CC events 
requires a precise measurement of all final state particles. 
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= 7.5 GeV 

E fl = 30 GeV 


Process 


L = 732 km 

L = 2900 km 




10 21 yr 

2.5 x 10 20 n~ 


0, CC 


39572 

35945 

Non-oscillated 

o* NO 


10775 

10688 

rates 

P e CC 


14419 

13911 


P e NC 


4501 

4830 

Oscillated 

P e P/x 

CC 

85 

42 

events (<5 = 7 t/2) 

Ox ^ "e 

CC 

248 

238 

Oscillated 

P e ^ P/x 

CC 

134 

72 

events (5 = 0) 

Ox V e 

CC 

370 

333 

Oscillated 

Pe ^ P/x 

CC 

134 

69 

events (5 = — vr/2) 

Ox ^e 

CC 

360 

323 


Table 1: Event rates for a 10 kton detector. The oscillation parameters are: Am 2 2 = 
3 x 10“ 3 eV 2 , Am 2 2 = 1 x 10 -4 eV 2 , sin 2 d 23 = 0.5, sin 2 d 12 = 0.5 and sin 2 26 i3 = 0.05. 





= 7.5 GeV 

= 30 GeV 


Process 


L = 732 km 

L = 2900 km 




10 21 //+ 

2.5 x 10 20 y+ 


P/x cc 


16472 

16054 

Non-oscillated 

P /t NC 


5319 

5636 

rates 

u e CC 


35279 

31399 


u e NC 


9268 

9245 

Oscillated 

^e Ox 

CC 

408 

389 

events (5 = n/2) 

P/x Pe 

CC 

123 

66 

Oscillated 

O ^ Ox 

CC 

403 

381 

events (5 = 0) 

P/X ^ Pe 

cc 

128 

71 

Oscillated 

^e Ox 

cc 

266 

273 

events (6 = —vr/2) 

P/x Pe 

cc 

84 

43 


Table 2: Same as Table 1, but n + decays. 
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e Vy Z/ e 


Right sign muons (p ): 



CC 



Vy ^ V T 

CC 

O' - 

-+ h') 

Ve 

NC 

O' - 

-»• p~) 

Vy 

NC 

O' - 

-»• h~) 


Wrong sign muons (/i + ): 


Ve Vy 

CC 



Ve V T 

CC 

(r+ - 


V e 

NC 

(7T + - 

-+P + ) 

Vy 

NC 

(7T + - 

-+P + ) 


Electrons: 


V e 

CC 

(pg -> e + ) 

Vy ^ V e 

CC 

(z/g -> e“) 

Vy V T 

CC 

0” -»■ e") 

Ve Vt 

CC 

(r + -> e + ) 


Wrong sign electrons (e ): 

Z/ e CC (z/ e * e ) X Pconf 

Vy CC (z/g * e ) X (1 pconf ) 

^ O- CC (r“ -> e“) x (1 - p conf ) 
_P e P T CC (t + —> e + ) x p conf 

Neutral Currents (no leptons): 

Uy NC (hadrons) 

z/ M v T CC (r“ —> hadrons) 

P e NC (hadrons) 

P e p T CC (r + — > hadrons) 


Table 3: Event classes in case of p beam. 


It is natural to classify the events in various classes depending on their final state 
configuration ||28|| . We illustrate them for the case of p~ stored in the ring. 

1. Right sign muons ( rsp ): the leading muon has the same charge as those 
circulating inside the ring. We include the non-oscillated Vy CC events, the Uy —> 
u T oscillations charged current events followed by r —> p and the background 
from hadron decays in neutral currents induced by all neutrino flavors. 

2. Wrong sign muons ( wsp ): the leading muon has opposite charge to those 
circulating inside the ring. Opposite-sign leading muons can only be produced 
by neutrino oscillations, since there is no component in the beam that could 
account for them. This includes P e —> Dy oscillations and P e —> P T oscillations 
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with r + —> p + decays. We also include hadron decays in neutral currents induced 
by all neutrino flavors. 

3. Electrons (e): events with a prompt electron ( both charges ) and no primary 

muon identified. Events with leading electron or positron are produced by the 
charged-current interactions of the following neutrinos: non-oscillated u e neu¬ 
trinos, Vp —> v e oscillations and v e —> D T or —» v T oscillations with t —> e 

decays 

4. Wrong sign electrons ( wse ): the leading electron has an identified charge 
which has the same sign as that of the muons circulating inside the ring. We 
define 

• the efficiency e e to identify the charge; 

• the probability p con f for charge confusion. 

Hence, the wrong sign electron sample is built from unoscillated P e ’s from the 
beam with a weight P(p e —> P e ) x e e xp con f and from muon neutrinos oscillations 
with the weight P(z/ M —> u e ) x e e x (1 - p con f )• 

5. No Lepton (Of): events corresponding to NC interactions or v T CC events 
followed by a hadronic decay of the tau lepton. Events with no leading electrons 
or muons will be used to study the v tl —> v T oscillations. These events can be 
produced in neutral current processes or in D e —> D T or z/ /t —> u T oscillations with 
r —> hadrons decays. 

Table |3| summarizes the different processes that contribute to the five event classes. 
Pconf is the electron charge confusion probability. 

For electron or muon charged current events, the visible event energy reconstructs 
the incoming neutrino energy and is therefore a very important variable to study the 
energy dependence of the oscillations. 

The last three classes can only be cleanly studied in a fine granularity detector. 

The energy loss of a minimum ionizing particle in liquid Argon is about 200 MeV/m, 
while the nuclear interaction length is 84 cm. That means that a muon with a mo¬ 
mentum of 0.5 GeV will travel about 2.5 meters before being absorbed, corresponding 
to about 3 nuclear interaction lengths. The probability that a pion will travel so long 
without interacting is small, and also detailed simulations show that at this level of 
background rejection pions do not represent a problem. 

Based on preliminary estimates^], we assume that for electrons with energy less 
than 7.5 GeV, the charge identification efficiency is 20% with a charge confusion prob¬ 
ability of 0.1%, unless otherwise noted. 
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11.2 Direct extraction of the oscillation probabilities 


From the visible energy distributions of the events, one can extract the oscillation 
probabilities. The visible energy of the events are plotted into histograms with 10 bins 
in energy. The v e —> v^ oscillation probability in each energy bin i can be computed as 


Ni(ws/i) — N®(wsp) 


[/i + decays) 


(59) 


where Nj{wsp) is the number of wrong-sign muon events in the i-th bin of energy, 
N^{wsy) are the background events in the i-th bin of energy, > p™*) is the ef¬ 
ficiency of the muon threshold cut in that bin, and iV°(e) is the number of electron 
events in the i-th bin of energy in absence of oscillations. The number of events corre¬ 
sponds to the statistics obtained from y + decays. A similar quantity for antineutrinos 
V l {v e — > P M ) will be computed with events coming from y~ decays. 

Similarly, the —> v e oscillation probability in each energy bin i can be computed 

as 





Nj{wse ) — N®(wse) 
e e (l — p C(m f)N^(rsn) 


[p decays] 


(60) 


where N^wse ) is the number of wrong-sign electron events in the i-th bin of energy, e e 
is the efficiency for charge discrimination, p con f the charge confusion, and Nf{rsp) is 
the number of right sign muon events in the i-th bin of energy in absence of oscillations. 
The number of events corresponds to the statistics obtained from p~ decays. A similar 
quantity for antineutrinos V t (P /t —> P e ) will be computed with events coming from p + 
decays. 

These binned probabilities could be combined in an actual experiment in order to 
perform direct searches of the effects induced by the 5-phase. 


11.3 Direct search for T-asymmetry 

For measurements involving the discrimination of the electron charge, we limit ourselves 
to the lowest energy and baseline configuration (A /t = 7.5 GeV and L = 732 km), since 
we expect the discrimination of the electron charge to be practically possible only at 
these lowest energies. 

The binned Ay(i) discriminant for neutrinos is defined as 

A T (i) = ViiVf. t —> Ve) - Vi{v e Vfj) (61) 

and a similar discrimant A r(i) can be computed for antineutrinos. 

These quantities are plotted for neutrinos and antineutrinos for three values of the 
5-phase (5 = +7t/2, 5 = 0 and 5 = — 7t/2) in Figure [L8]. The errors are statistical and 
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correspond to a normalization of 10 21 muon decays and a baseline of L — 732 km. A 
20% electron efficiency with a charge confusion probability of 0.1% has been assumed. 
The full curve corresponds to the theoretical probability difference. 

A nice feature of these measurements is the change of sign of the effect with re¬ 
spect of the change 5 —> — 5 and also with respect to the substitution of neutrinos by 
antineutrinos. These changes of sign are clearly visible and would provide a direct, 
model-independent, proof for T-violation in neutrino oscillations. 

In order to cross-check the matter behavior, one can also contemplate the CPT- 
discriminants defined as 


A cPT^i) = Viivfj, —> u e ) - Vi{v e —> uf) 

A cptH) = Vi{y e -> Vp) - Viiy^ -> V e ) (62) 

These quantities are plotted in Figure [19], with the same assumptions as in Figure [L8|. 
The full curve corresponds to the theoretical probability difference. As expected, these 
quantities are independent from the 5-phase and probe only the matter effects. The 
change of sign of the effect with respect to the substitution of neutrinos by antineutrinos 
is clearly visible. 

It should be however noted that in the case of the CPT discriminant, the statisti¬ 
cal power is rather low, since this measurement combines the appearance of electrons 
(driven by the efficiency for detecting the electron charge) and involves antineutri¬ 
nos, which are suppressed by matter effects. Hence, the statistical power is reduced 
compared to the T-discriminant. 

11.4 Direct search for CP-asymmetry 

In the direct search for the CP-asymmetry, we rely only on the appearance of wrong- 
sign muons. We compare in this case the two energy and baselines options. 

The binned A cpif) discriminant for the shortest baseline L = 732 km, E^ = 
7.5 GeV and longest baseline L = 2900 km, Ef = 30 GeV (lower plots) for three values 
of the h-phase (5 = +7t/2, 5 = 0 and 5 = — vr/2) are shown in Figure The errors are 
statistical and correspond to a normalization of 10 21 (2.5 x 10 20 ) for L = 732(2900) km. 
The full curve corresponds to the theoretical probability difference. The dotted curve 
is the theoretical curve for 5 = 0 and represents the effect of propagation in matter. 

As was already pointed out, the A C p does not vanish even in the case <5 = 0, 
since matter introduces an asymmetry. At the shortest baseline (L = 732 km), these 
effects are rather small, as illustrated in Figure 0 This has the advantage that the 
observed asymmetry would be positive for 5 > 0, but would still change sign in the 
case 5 ~ — n/2. In the fortunate case in which Nature has chosen such a value for 
the 5-phase, the observation of the negative asymmetry would be a striking sign for 
CP-violation, since matter could never produce such an effect. 
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For other values of the 5-phase, the effect is positive. It is also always positive at 
the largest baseline L = 2900 km, since at those distances the effect induced by the 
5-phase is smaller than the asymmetry introduced by the matter. 


11.5 Comparison of two methods 


The binned A r(i) and A cp(i) discriminant can be used to calculate the y 2 significance 
of the effect, given the statistical error on each bin. We compute the following y 2 ’s: 


Xt ~ 


(Ar(t, 5) — A r(i, 5 = 0)) 2 (A r(i, 5) — A r(i, 5 — 0)) 


<7(A r (i,5))= 


+ 


a(A T (i,S)y 


(63) 


where a(A T (i,6)) is the statistical error in the bin. Since A r (5) and A T (S) use inde¬ 
pendent sets of data, the global \t can be computed in this way, as linear sum of both 
contributions. 

Similarly, the y 2 of the CP-asymmetry is 


X 2 cp = 

i 


(A cp(i, $) — Acp{i, 5 — 0)) 2 
o-(Ac- P (i,5)) 2 


(64) 


We study the significance of the effect as a function of the solar mass difference 
Ara^i, since the effect associated to the 5-phase will decrease with decreasing A 
values. We consider the range compatible with solar neutrino experiments, 10 -5 < 
Am 2 i < 10- 4 eV 2 . 

The exclusion regions obtained at the 90%C.L. (defined as Ay 2 = +1.96) in the 
5-phase vs Am^ plane are shown in Figure 0 The rest of parameters are fixed to 
the reference values (Am| 2 = 3 x 1CT 3 eV 2 , Am 21 = 1 x 10“ 4 eV 2 , sin 2 023 = 0.5, 
sin 2 0 12 = 0.5, and sin 2 20 13 = 0.05). An electron charge confusion probability of 0.1% 
and an electron detection efficiency of 20% has been assumed. The normalizations 
assumed are 10 21 and 5 x 10 21 muon decays with energy A /t = 7.5 GeV and a baseline 
of L — 732 km. 


The results are very encouraging. With 10 21 muon decays, the region Am 21 > 
4 x 10~ 5 eV 2 is covered. For 5 x 10 21 muons, this region extends down to 2 x 10 -5 eV 2 . If 
we consider that the value of Am 2 ] is known and that it has a value of Am 21 = 10~ 4 eV 2 , 
one can constrain the values of the 5-phase within the range |5| < 0.57 or |5| > 2.6 for 
10 21 muons and |5| < 0.12 and |5| > 3.0 for 5 x 10 21 muon decays at the 90%C.L. 


Figure 22 shows the same exclusion plot but comparing the two considered L/E v 
values. For L = 730 km, the normalizations are 10 21 and 5 x 10 21 muons, while for 
L = 2900 km the considered fluxes are four times smaller. Contrary to what happens 
at short baselines, where a nice symmetry between +5 and —5 is observed, at L = 
2900 km matter effects introduce a clear asymmetry between the two excluded regions. 
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We conclude that an exhaustive direct, model-independent exploration of the 5- 
phase, within the full range 10 -5 < Am 21 < 10~ 3 4 5 eV 2 requires an intensity of 5 x 10 21 
muon decays of each sign. 

As explained in Section |11.2| , the value of the —> v e oscillation probability (and 
hence the Xt discriminant) depends on both the electron detection efficiency, e e , and 
the electron charge confusion, p con f (see equation |60|). In Figure ^3] we show the iso¬ 
curves at la, 2a and 3er levels of significance^ that can be obtained for S = +7t/ 2 
and —7 t/ 2. For instance, in the case of (e e , p con f)=( 20%,1%) one could exclude both 
5 = ±7t/ 2 at 2a level. Looking at the shape of the iso curves we can conclude that, 
even if the result depends on both parameters, it is much more sensitive to the value 
of the charge confusion, which should be smaller than ~0.1% for efficiencies of ~20% 
to have an exclusion power above 3 a. 

11.6 Comparison with the fit of the visible energy distribu¬ 
tions 

The most effective way to fit the oscillation parameters is to study the visible energy 
distribution of the four classes of events^, since assuming the unoscillated spectra are 
known, they contain direct information on the oscillation probabilities. 

Of course, for electron or muon charged current events, the visible energy recon¬ 
structs the incoming neutrino energy. In the case of neutral currents or the charged 
current of tau neutrinos, the visible energy is less than the visible energy because of 
undetected neutrinos in the final state. The information is in this case degraded but 
can still be used. 

In case no CP-violation is observed, the result of the fit in terms of 2-dimensional 
90% C.L. contours in the Am 2 2 — 6 plane is shown in Figure [24|. The Y axis spans the 
A'/7i 2 2 range allowed by LMA solar neutrinos, and the whole range — n < 5 < tt has 
been consider. For each pair of values (Am 2 2 ,5), the fit to the reference distribution^ 
was performed fixing all parameters (lower curves) and leaving # 13 free (upper curves). 
The result for the two considered baselines and energies are shown. 

As expected, this method is clearly more powerful than the direct search for CP- 
asynnnetry discussed in Section |11.4| . The comparison between both results reveals 
that the non-excluded zones near the extremes (|<5| ~ tt) in Figure [H] are now covered 
and can be explored with this second method (Figure |24|). However, even if higher 
exclusive, the fit of the visible energy distributions requires a good knowledge of the 
oscillation parameters, while the direct search for C'P-asynimetry is essentially model 

3 The ncr region is reached when ^J\at + — n • 

4 In this case, we do not use the information coming from the electron charge. 

5 The reference values are Arn| 2 = 3 x 10 -3 eV 2 , sin 2 d 2 3 = 0.5, sin 2 0i 2 = 0.5, sin 2 20 i 3 = 0.05 and 
<5 = 0 . 
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independent. Nevertheless, even assuming that the value of On is not known precisely 
(upper curves in Figure |4|) the excluded regions extend beyond the limits set by the 
first method for |h| ~ 7r. 

The result obtained with the full simulation fully supports all the above phenomeno¬ 
logical considerations: for baselines such that the maximum of the interesting effect lies 
well below the MSW resonance, the destructive effect of matter almost does not play 
any role, and for a given machine power, if L/E^ is kept constant, the same sensitivity 
is reached. 


12 Summary and Conclusions 

In this document, we have discussed general strategies to detect CP-violation effects 
related to the complex phase 5 of the neutrino mixing matrix, in the framework of a 
neutrino factory. Here is the summary and conclusions: 

• In order to directly compare effects at different energies and baselines, we have 
defined the “rescaled probability”: P[v a -w vp]E v ,L) x -jj. It approximately weighs 
the probability by the neutrino energy spectrum of the neutrino factory and by the 
attenuation of the neutrino flux with the distance (Section [!]). 

• In vacuum, the region of the “first maximum” on the —> v e oscillation proba¬ 
bility yields E™ ax ~ 2 GeV at 730 km, ~ 8 GeV at 2900 km and ~ 20 GeV at 7400 km 
for Am§ 2 = 3 x 10~ 3 eV 2 . For neutrinos propagating through matter, the oscillation 
probabilities are distorted and the resonant energy is E xes ~ 14.1, 12.3 and 10.7 GeV 
for Am| 2 = 3 x 10~ 3 eV 2 and p equal to 2.7, 3.2 and 3.7 g/cm 3 respectively. 

• The most favorable choice of neutrino energy E v and baseline L is in the region of 
the “first maximum” given by {L/E v ) max ~ 400 km/GeV for |A/n 2 2 | = 3 x 10~ 3 eV 2 . 
As discussed in Section |64], if one wants to study oscillations in this region one has to 
require that the energy of the “first-maximum” be smaller than the MSW resonance 
energy: 2\/2GFn e E™ ax < Am§ 2 cos 2# 13 . This fixes a limit on the baseline L max «A000 
km beyond which matter effects spoil the sensitivity. 

This implies that we concentrate on medium L/E v . 

• We have considered three quantities which are good discriminators for a non¬ 
vanishing h-phase (Section |6|): A $, A cp and At- 

A s can provide excellent determination of the phase, limited only by statistics, but 
needs a precise knowledge of the other oscillation parameters and possible correlations 
with the d 13 at high energy can appear. 

A cp has the advantage that involves the appearance of wrong-sign muons only, 
experimentally easy to detect, but the disadvantage of involving both, neutrinos and 
antineutrinos, which requires a good understanding of the effects related to matter. 
Matter effects will largely spoil A cp at long distances. 
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A t is not affected by matter effects because it involves only neutrinos, but it requires 
to discriminate the electron charge, which is experimentally very challenging. A non¬ 
vanishing Ax would be a direct proof for a non-vanishing J-phase. 

We have demonstrated that the discriminants of the 5-phase in vacuum scale with 
L/E u (in matter, this scaling remains valid for L smaller than ~5000 km). 

• In one hand, because of the linear rise of the neutrino cross-section with neutrino 
energy at high energy, the statistical significance of the effect scales with E v , so it 
grows linearly with L (for L/E u constant). 

On the other, the study of A^ requires the electron charge identification. For a 
constant (L/E u ) ratio, the need of low energy electrons points towards lower-energy 
beams and shorter distances. 

• As an example, we have considered two concrete cases keeping the same L/E M 
ratio: (L = 730 km, E^ = 7.5 GeV) and (L = 2900 km, Ef = 30 GeV). We have 
classified the events in five classes (right and wrong sign muons, right and wrong sign 
electrons and no leptons) assuming a 10 kT Liquid Argon detector with an electron 
charge confusion of 0.1% . We have computed the exclusion regions in the A m\ 2 — 
5 plane using the A cp and A t discriminants, obtaining a similar excluded region 
provided that the electron detection efficiency is ~20%. The A m\ 2 compatible with 
the LMA solar data can be tested with a flux of 5xl0 21 muons. 

• Finally, we have computed the exclusion regions in the A m\ 2 — 5 plane fitting 
the visible energy distributions. This method, more powerful than the previous one 
but not model independent, extends the excluded regions up to values of |<5| close to 
7r, even when is left free. 


References 


[1] Y. Fukuda et al. (Kamiokande Collaboration), Phys. Lett. B 335, 237 (1994); R. 
Becker-Szendy et al. (IMB Collaboration), Nucl. Phys. B (Proc. Suppl.) 38, 331 
(1995); S. Fukuda et al. (Super-Kamiokande Collaboration), |hep-ex/0009001| ; T. 
Kajita (Super-Kamiokande Collaboration), Talk presented at NOW2000, Otranto, 
Italy, September 2000 ([http:/ /www.ba.infn.it/~now2000|) ; W.W.M. Allison et al. 
(Soudan 2 Collaboration), Phys. Lett. B 449, 137 (1999); F. Ronga (MACRO 
Collaboration), Talk presented at NOW2000, Otranto, Italy, September 2000 
(|http:/ /www. ba.infn.it/ ~now2000|). 


[2] B. Pontecorvo, J. Expt. Theor. Phys. 33, 549 (1957) [Sov. Phys. JETP 6, 429 
(1958)]; B. Pontecorvo, J. Expt. Theor. Phys. 34, 247 (1958) [Sov. Phys. JETP 
7, 172 (1958)]; Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28 
(1962) 870; B. Pontecorvo, J. Expt. Theor. Phys 53 (1967) 1717; V. Gribov and 
B. Pontecorvo, Phys. Lett. B 28, 493 (1969). 


31 









[3] 

[4] 

[5] 

[ 6 ] 


[7] 

[ 8 ] 

[9] 

[ 10 ] 


[ 11 ] 

[ 12 ] 

[13] 


“E362 Proposal for a long baseline neutrino oscillation experiment, using KEK-PS 
and Super-Kamiokande”, February 1995. 

H. W. Sobel, proceedings of Eighth International Workshop on Neutrino Tele¬ 
scopes, Venice 1999, vol 1 pg 351. 


E. Abies et al. [MINOS Collaboration], “P-875: A Long baseline neutrino oscilla¬ 
tion experiment at Fermilab,” FERMILAB-PROPOSAL-P-875. 

The MINOS detectors Technical Design Report, NuMI-L-337, October 1998. 


K. Kodama et al. [OPERA Collaboration], “OPERA: a long baseline v T appearance 
experiment in the CNGS beam from CERN to Gran Sasso”, CERN/SPSC 99-20 
SPSC/M635 LNGS-LOI 19/99. 

ICARUS Collab., ’ICARUS-II. A Second-Generation Proton Decay Experiment 
and Neutrino Observatory at the Gran Sasso Laboratory Proposal Vol. I & II, 
LNGS-94/99, 1994. ICARUS Collab. “A first 600 ton ICARUS Detector Installed at 
the Gran Sasso Laboratory”, Addendum to proposal, LNGS-95/10 (1995). F. Ar- 
neodo et al. [ICARUS Collaboration], “Cloning of T600 modules to reach the design 
sensitive mass”, LNGS-EXP 13/89 add. 2/01, Nov. 2001 and references therein. 
Available at | http://www. cern. ch/icaru 4 


F. Arneodo et al. [ICARUS collaboration], “The ICARUS experiment, a second- 
generation proton decay experiment and neutrino observatory at the Gran Sasso 
laboratory,” |hcp-cx/0103008| . 


F. Arneodo et al. [ICARUS and NOE Collaborations], “ICANOE: Imaging 
and calorimetric neutrino oscillation experiment: a proposal for a CERN-GS 
long baseline and atmospheric neutrino oscillation experiment.” LNGS-P21/99; 
CERN/SPSC 99-25; SPSC/P314, September 99. 


S. Geer, Phys. Rev. D 57, 1998 (6989) 


A. de Rujula, M. B. Gavela and P. Hernandez, Nucl. Phys. B 547, 1999 (21); 
V. Barger, S. Geer, R. Raja and K. Whisnant, Phys. Rev. D62, 013004 (2000) 
hcp-ph/991 152411 ; Phys. Rev. D62, 073002 (2000) [|icp-ph/0003184|l ; A. Bueno, 


M. Campanclli and A. Rubbia, Nucl. Phys. B573 (2000) 27; V. Barger, S. Geer 
and K. Whisnant, Phys. Rev. D 61, 2000 (053004); M. Freund, M. Lindner, S. 
T. Petcov and A. Romanimo, Nucl. Phys. B578, 27 (2000) | [hcp-ph/991 245 7|| ; A. 
Cervera et al. , Nucl. Phys. B579, 17 (2000) |[hcp-ph/0002108|| . 


A. Bueno, M. Campanclli and A. Rubbia, |hcp-ph/9809252| ; CERN-EP/98-140. 


M. Apollonio et a/.[CHOOZ Collaboration], Phys. Lett. B 466, 1999 (415) Phys. 
Lett. B 420, 1998 (397) 


A. Cervera et al ., Nucl. Phys. B579, 17 (2000) |[licp-ph/0002108 


32 


























[14] Information on the neutrino factory studies and mu collider collaboration at BNL 
can be found at | http://www.cap.bnl.gov/mumu/ 

Information on the neutrino factory studies at FNAL can be found at 
http://www.fnal.gov/projects/muon_collider/ 


Information on the neutrino factory studies 
| http://muonstoragerings. cern. ch/Welcome, html/ 


at CERN can be found at 


[15] A. Cervera, A. Donini, M. B. Gavela, J. J. Gomez Cadenas, P. Hernandez, O. Mena 
and S. Rigolin, Nucl. Phys. B 579, 17 (2000) ||hcp-ph/0002108| . 

[16] P. Lipari, “CP violation effects and high energy neutrinos”, |hcp-ph/0102046| . 


[17] A. Bueno, M. Campanelli and A. Rubbia, Nucl. Phys. B 589, 577 (2000) |[hcp- 
| ph/0005007[ . 


[18] T. Ota and J. Sato, “Matter profile effect in neutrino factory”, |hcp-ph/0011234 . 

[19] M. Koike and J. Sato, Phys. Rev. D 61, 073012 (2000) ||hcp-ph/9909469| . 


[20] J. Arafune, M. Koike and J. Sato, Phys. Rev. D 56, 3093 (1997) ||hcp-ph/9703351 
J. Arafune and J. Sato, Phys. Rev. D 55, 1653 (1997) 


[21] J. Sato, “CP and T violation in (long) long baseline neutrino oscillation experi¬ 
ments” , |hcp-ph/9910442| . 

[22] H. W. Zaglauer and K. H. Schwarzer, Z. Phys. C 40, 1988 (273) 


[23] V. Barger, K. Whisnant, S. Pakvasa and R. J. N. Phillips, Phys. Rev. D 22, 1980 
(2718) 

[24] P. Lipari, Phys. Rev. D 61, 113004 (2000) ||hcp-ph/9903481| . 

[25] P. F. Harrison and W. G. Scott, Phys. Lett. B 476, 349 (2000) ||hcp-ph/9912435| ■ 

[26] A. Romanino, Nucl.Phys. B 574 675-690 (2000) ||hcp-ph/9909~425[ 

[27] A. Rubbia et al. , talk given at NUINT01, KEK (Japan), Dec. 2001. 


[28] A. Bueno, M. Campanelli and A. Rubbia, |hcp-ph/9808485| ETHZ-IPP-98-05. 


[29] S. Navas et al, “Prospects for CP and T violation measurements with a mag¬ 
netized Liquid Argon detector”, talk given at the Workshop on a Massive Un¬ 
derground Neutrino Detector with Leptonic Charge-Sign Discrimination, Sitges 
(Spain), Nov. 2001. A. Rubbia, “Neutrino factories: Detector concepts for stud¬ 
ies of CP and T violation effects in neutrino oscillations,” |arXiv:hcp-ph/0106088| , 
talk given at 9th International Symposium on Neutrino Telescopes, Venice (Italy), 
March 2001. 


33 
















































L=730 km 



1 10 io 2 

E v (GeV) 

L = 2900 km 



1 10 io 2 

E v (GeV) 


Figure 1: Probability for u e — > oscillations in vacuum for two baselines L = 730 km 
and 2900 km as a function of neutrino energy E v . The probabilities are computed for 
three values of the 5-phase: 5 = 0 (line), 5 = + 7 t /2 (dashed), 5 = — 7 t /2 (dotted). 
The other oscillation parameters are Am | 2 = 3 x 10 “ 3 eV 2 , A= 1 x 10 “ 4 eV 2 , 
sin 2 6*23 = 0.5, sin 2 0 12 = 0.5, and sin 2 20 13 = 0.05. 
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Figure 2: Rescaled probability (see text) for u e 
baselines L = 730 km (line), 2900 km (dashed) and 7400 km (dotted) as a function of 
neutrino energy E v . The oscillation parameters are A m\ 2 = 3 x 10 
1 x 10 -4 eV 2 , sin 2 d 2 3 = 0.5, sin 2 0 12 = 0.5, sin 2 2d 13 = 0.05 and <5 = 0. 
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Figure 3: Rescaled probability (see text) for u e —> v lx oscillations in vacuum for three 
baselines L = 730 km (line), 2900 km (dashed) and 7400 km (dotted) as a function 
of neutrino energy E u . For each baseline, the probabilities are computed for three 
values of the 5-phase: 5 = 0 (upper), 5 = +7 t/2 (middle), 5 = —7 t/2 (lower). The 
oscillation parameters are Am | 2 = 3 x 10 ~ 3 eV 2 , Am 2 ] = 1 x 10 -4 eV 2 , sin 2 0 2 3 = 0.5, 
sin 2 9i2 = 0.5, sin 2 2d 13 = 0.05. 
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Figure 4: Oscillation probability for v e —> oscillations for a baseline L = 730 km 
as a function of neutrino energy E v . The probabilities are computed for neutrinos 
in matter (full line) and in vacuum (dotted line) and for three values of the h-phase: 
8 — 0, 8 — + 7 t /2 and 8 = — 7 t/ 2 . The oscillation parameters are A m\ 2 = 3 x 10 ~ 3 eV 2 , 
A= 1 x 10 ” 4 eV 2 , sin 2 6*23 = 0.5, sin 2 d 12 = 0.5, sin 2 26* 13 = 0.05. 
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Figure 5: Same as Figure [| but for a baseline L = 2900 km. 
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Figure 6 : Rescaled probability (see text) for v e —> oscillations for a baseline 
L = 730 km as a function of neutrino energy E v . The probabilities are computed 
for neutrinos in matter, in vacuum (dotted line, same for neutrinos and antineutrinos) 
and for antineutrinos in matter. The oscillation parameters are = 3 x 10 ~ 3 eV 2 , 

Amlx = 1 x 10 ” 4 eV 2 , sin 2 6*23 = 0.5, sin 2 0 12 = 0.5, sin 2 26*13 = 0.05 and 5 = 0. 


39 







-6 


x 10 


I >-2900 km, p=3.2 g/cm 3 



Figure 7: Same as Figure |6] but for a baseline L = 2900 km. 
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Figure 8: Same as Figure |6] but for a baseline L = 7400 km. 
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Figure 9: Same as plots in Figure |S| but only for neutrinos and as a function of L/E u 
for three baselines L = 730 km, 2900 km, 7400 km and in vacuum (dashed line, 
independent of baseline). The oscillation parameters are A m\ 2 = 3 x 10 -3 eV 2 , A= 
1 x 10 -4 eV 2 , sin 2 6b 3 = 0.5, sin 2 0 12 = 0.5, sin 2 20 13 = 0.05 and 5 = 0. 
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Figure 10: The rescaled As discriminant (see text for definition) as a function of the 
L/E v ratio, computed for neutrinos propagating in matter at three different baselines 
L = 730 km, 2900 km and 7400 km (full lines), and for neutrinos propagating in 
vacuum (dashed line). 
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Figure 11: Rescaled probability (see text) for v e —*■ oscillations for a baseline 
L = 730 km as a function of neutrino energy E v . The probabilities are computed 
for neutrinos in matter (full line) and in vacuum (dotted line), and for three val¬ 
ues of the 5-phase: S = 0, S = +7t/2, 5 = — tt/2. The other oscillation parameters 
are A m \ 2 = 3 x 10 -3 eV 2 , Am| x = 1 x 10~ 4 eV 2 , sin 2 0 2 3 = 0.5, sin 2 612 = 0.5, 
sin 2 26*13 = 0.05. 
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Figure 12: Same as Figure 11 but for a baseline L 


2900 km. 
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Figure 13: The rescaled A t discriminant (see text for definition) as a function of the 
L/E v ratio, computed for neutrinos propagating in matter at three different baselines 
L = 730 km, 2900 km and 7400 km, and also for propagation vacuum (independent of 
baseline). Two sets of curves are represented, corresponding to 5 = +7 t/2 (upper ones) 
and 6 = —7 t/2 (lower ones). The other oscillation parameters are Am 2 2 = 3 x 10 ~ 3 eV 2 , 
Am 2 ! = 1 x 10 -4 eV 2 , sin 2 6*23 = 0.5, sin 2 9 12 = 0.5 and sin 2 26* 13 = 0.05. 
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Figure 14: The rescaled A^ discriminant (see text for definition) as a function of 
the L and E v , computed for neutrinos propagating in vacuum (upper) and in matter 
(lower) for S = tt/2. The other oscillation parameters are Am § 2 = 3 x 10 -3 eV 2 , 
Amli = 1 x 1CT 4 eV 2 , sin 2 0 2 3 = 0.5, sin 2 0\ 2 = 0.5 and sin 2 26*i 3 = 0.05. 
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Figure 15: The rescaled A cp discriminant (see text for definition) as a function of the 
neutrino energy E u , computed for neutrinos propagating in matter at three different 
baselines L = 730 km, 2900 km and 7400 km. Three sets of curves are represented, 
corresponding to 6 = +7 t/2 (thick lines), 5 = —7 t/2 (thin lines) and 5 = 0. The other 
oscillation parameters are A m\ 2 = 3 x 10 ~ 3 eV 2 , Am)] = 1 x 10 -4 eV 2 , sin 2 623 , = 0.5, 
sin 2 9 12 = 0.5 and sin 2 29 13 = 0.05. 
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Figure 16: Same as Figure [l5]but as a function of the L/E v ratio. In addition to the 
three different baselines, the case of neutrinos propagating in vacuum (independent of 
baseline) is also shown for comparison. 
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Charge confusion = 0.1 % , Electron efficiency = 20 % 
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Figure 17: Visible energy distribution for wrong-sign muons (left) and wrong-sign 
electrons (right) normalized to 10 21 muon decays. The electron efficiency e e is assumed 
to be 20% and charge confusion probability is set to 0.1%. Three sets of curves are 
represented, corresponding to S = +7t/ 2 (dashed line), 5 = 0 (full line) and S = —7t/ 2 
(dotted line). The background contribution from tail decays is also shown. The other 
oscillation parameters are A m\ 2 = 3 x 10 -3 eV 2 , Am|i = 1 x 10 -4 eV 2 , sin 2 0 2 3 = 0.5, 
sin 2 9 12 = 0.5 and sin 2 26 *i 3 = 0.05. 
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Figure 18: Direct T-violation: Binned A r(i) discriminant for neutrinos (upper plots) 
and antineutrinos (lower plots) for three values of the h-phase: S = +7t/ 2, <5 = 0 and 
5 = —7t/ 2. The errors are statistical and correspond to a normalization of 10 21 muon 
decays and a baseline of L = 732 km. A 20% electron efficiency with a charge confusion 
probability of 0.1% has been assumed. The change of sign of the effect with respect of 
the change 6 —> —6 and the substitution of neutrinos by antineutrinos is clearly visible. 
The full curve corresponds to the theoretical probability difference. 
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Figure 19: Direct CPT-violation: Binned A cpt{}) (upper plots) and Acpt(}) (lower 
plots) for three values of the 5-phase: 5 = +7t/2, 5 = 0 and 5 = — 7t/2. The errors 
are statistical and correspond to a normalization of 10 21 muon decays and a baseline 
of L = 732 km. A 20% electron efficiency with a charge confusion probability of 0.1% 
has been assumed. The change of sign of the effect with respect to the substitution 
of neutrinos by antineutrinos is clearly visible. The full curve corresponds to the 
theoretical probability difference. 
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Figure 20: Direct CP-violation: Binned A cp(i) discriminant for the shortest baseline 
L = 732 km, = 7.5 GeV (upper plots) and longest baseline L = 2900 km, P M = 
30 GeV (lower plots) for three values of the 5-phase: 5 = +7t/2, 5 = 0 and 5 = — 7t/2. 
The errors are statistical and correspond to a normalization of 10 21 (2.5 x 10 20 ) for L = 
732(2900) km. The full curve corresponds to the theoretical probability difference. The 
dotted curve is the theoretical curve for 5 = 0 and represents the effect of propagation 
in matter. 
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Figure 21: Exclusion region at the 90%C.L. (x 2 > / \mm + 1-96) in the (5-phase vs Am|i 
plane. Two regions obtained with the A cp and the sum of A t and A t discriminants 
are shown. A charge confusion probability of 0.1% and an electron efficiency of 20% 
has been assumed. The result is shown for two normalizations, 10 21 and 5 x 10 21 muon 
decays of each type with energy E M = 7.5 GeV and a baseline of L = 732 km. The 
reference oscillation parameters are Am 2 2 = 3 x 10 -3 eV 2 , sin 2 6b 3 = 0.5, sin 2 0 12 = 0.5 
and sin 2 2ffi 3 = 0.05. 
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Figure 22: Exclusion region at the 90%C.L. (y 2 > Xmin + 1-96) in the 5-phase vs 
A rri^i from the A cp discriminant. The result is shown for two baselines and two 
normalizations: (L = 732 km, £0 = 7.5 GeV, </>o = 10 21 /i flux) and (L = 2900 km, 
E M = 30 GeV, 0o — 2.5 x lO 20 ^ flux). The reference oscillation parameters are Am § 2 = 
3 x 10 -3 eV 2 , sin 2 6*23 = 0.5, sin 2 #i 2 = 0.5 and sin 2 20 13 = 0.05. 
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Figure 23: The la, 2a and 3a y 2 contours of the combined A^ and A^ discriminants 
as a function of the electron efficiency and the probability for charge confusion. The 
result is shown for 5 = +7t/2 and 6 = — 7t/2. The other oscillation parameters are 
A 777.32 = 3 x 10- 3 eV 2 , Am|i = 1 x 1CT 4 eV 2 , sin 2 0 2 3 = 0.5, sin 2 0 i2 = 0.5 and 
sin 2 2013 = 0.05. 
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Figure 24: 90% C.L. sensitivity on the CP-phase S as a function of f° r the two 

considered baselines. The reference oscillation parameters are Am 2 2 = 3 x 10 ~ 3 eV 2 , 
sin 2 $23 = 0.5, sin 2 612 = 0.5, sin 2 20 13 = 0.05 and S = 0. The lower curves are made 
fixing all parameters to the reference values while for the upper curves #13 is free. 
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